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Abstract. We study the instability of a moving spinor Bose-Einstein 
condensate when the speed of flow reaches the critical velocity. This 
we identify on the basis of Landau's criterion, i.e. the velocity above 
which some elementary excitation energy becomes negative. We show 
that the first-to-become unstable excitations are spin-carrying quasi- 
particles. We also discuss the temperature dependence of the critical 
velocity in a more advanced mean-field approximation. 



1 Introduction 



The nature of superfluidity of Bose systems is the result of a remarkable interplay 
between Bose-Einstein condensation and particle-particle interaction. The complex 
phenomena of superfluids had been explained by a remarkably simple two fluid model 
[T]. One of the key concepts is the critical velocity, that is, the limit velocity of an 
obstacle immersed in the superfluid and moving against it without exciting density 
waves, therefore the motion is exempt from damping. In Landau's original argumen- 
tation for the critical velocity only Galilean invariance and energy and momentum 
conservation were taken into account besides a specific dispersion relation of the 
elementary excitations [I]. Later it was shown by Bogoliubov that such dispersion 
relation is indeed the consequence of Bose-Einstein condensation and the critical ve- 
locity is the speed of sound of density waves. In gases of ultracold atoms the critical 
velocity has been measured experimentally both for bosons [2 and also for fermions 
[3j. Though some doubt due to the inhomogeneity of the condensate and the finite 
size of the disturbing potential remains, the experiments confirm a speed of sound 
around that of the Landau criterion. Recently Navez and Graham showed for a spin- 
less Bose gas in a dynamically consistent Hartree-Fock-RPA approximation [4 that 
at the critical velocity not only the quasiparticle energy becomes zero but also the 
damping rate changes sign signaling the onset of a dynamical instability. 

Here we show for antiferromagnetic spin-1 Bose gases that the critical velocity 
is the speed of spin-wave instead of the sound wave. In typical experiments with 
ultracold gases the spin waves are usually two orders of magnitude slower than sound 
waves. Therefore the critical velocity is much smaller than in a scalar gas. A direct 
physical consequence is that when a spinor Bose condensate moves faster than the 
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speed of the spin wave, it emits a cone of spin polarized excitations instead of density 
waves, which can lead to spin texture formation and can be measured in experiments. 



2 The formulation of the problem of an antiferromagnetic spin-1 
Bose gas with a moving condensate 

We consider an ultracold, dilute, homogeneous, spin-1 Bose gas in a homogeneous 
magnetic field. The interparticle interaction is modeled by s-wave scattering, i.e. we 
neglect the relatively weak dipolar interaction of the gas. The grand-canonical Hamil- 
tonian of the system reads as 



+ W ^ aUki)at(k2)l/;/'a.(k3)as'(k4), (1) 

ki+k2=k3+k4 

where (k) is the annihilation operator of plane wave states with momentum k and 
spin projection r. The spin index r refers to the eigenvalue of the z-component of the 
spin operator and can take values from 0, — . Correspondingly = diag(l, 0, —1) 
is a 3x3 diagonal matrix. In Eq. ([T]) Ck = h?k'^ / {2M) refers to the kinetic energy of 
an atom, /i to the chemical potential, g to the gyromagnetic ratio, /xb to the Bohr 
magneton, B to the modulus of the homogeneous magnetic field. V is the volume of 
the system and V^j," the amplitude of the two particle interaction, given for spin-1 
bosons by [51617) : 

V;'/ = cJrsSr's' + Cs{¥)rs{'F)r's'. (2) 

with c„ — A'nh^{aQ + 2a2)/(3Af) , and Cg — 47rfi^(a2 — ao)/(3M). The parameters ao 
and a2 are the scattering lengths in the total hyperfine spin channel zero and two, 
respectively. The constant c„ > 0, while Cg can both be positive or negative, depending 
on the relative values of oq and 02. For > a zero net spin is energetically favorable 
(in the absence of a magnetic field). For this reason systems with Cs > are called 
antiferromagnetic or polar gases [5j. For instance the ultracold gas of ^"^Na in the 
f — I hyperfine state is antiferromagnetic [8]. 

Even in the presence of a homogeneous magnetic field the z component of the 
total spin is a conserved quantity. This conservation law can be resolved similarly as 
that for the particle number [i.e., with the introduction of a Lagrange multiplier in 
the Hamiltonian ([T])] . This multiplier shows up in the same way as the magnetic field 
does, therefore an effective magnetic field can be introduced as a sum of the external 
magnetic field plus the Lagrange multiplier. In the following, B will mean this kind 
of effective field even when the physical magnetic field is zero. In the following, we 
assume that either such is the case, or the external magnetic field is so small that the 
quadratic Zeeman effect can be neglected. 

In the presence of a Bose- Einstein condensate moving with velocity vq — ftko/M 
the annihilation and creation operators have a nonzero mean value ar.(ko) = ^/NcCr, 
where Nc is the number of atoms in the condensate, and Cr = (C-i-)Co;C-) is the 
condensate spinor normalized to unity 5 . It is convenient to define new annihilation 
operators with the canonical transformation br(k) = ar(k) — •\/-/VcCr'^k.ko which have 
zero mean values for all k. The finite temperature Green's function of the system is 
introduced as 

^-(k,r) = -(rA"(k,r)6f(k,0)), (3) 
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where r is the imaginary time, T,- is the r ordering operator and the Greek indices take 
the values ±1 and are introduced for a shorthand notation, namely b^(k, r) = t) 
and 6^^(k, r) — 6j!(2ko — k, — r). The Green's function can be conveniently arranged 



into a matrix form 
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(4) 



In Matsubara representation the inverse of the Green's function is expressed from the 
Dyson-Beliaev equations with the help of the free propagator and the self-energy 



G ^(k,ia;„) = Gq ^(k, itj„) - Z(k, ia;„). 



(5) 



The self-energy matrix Z is arranged from its components in a similar way than the 
Green's function Q. The inverse of the free propagator is given by the diagonal 
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(6) 



with E = iLOn — fi~^(ek — /i), E = —iuin — fi~^(e2ko-k — and ujl = hr^g^iBB the 
Larmor frequency related to the effective magnetic field. The self-energy Z(k, iuJn) is 
to be chosen according to some approximation scheme. The dispersion relation w(k) 
of the quasiparticles is provided by the poles of the retarded Green's function, that is 



dct G-'^ {k, oj + ir]) = 0. 



(7) 



To obtain the retarded Green's function we analytically continue in the frequency 
variable iw„ — > w + iry (with 77 — > +0) through the upper half complex plane. In 
the following sections we are going to evaluate the self-energy first in the Bogoliubov 
approximation valid at zero temperature, and later in the Hartree-RPA approximation 
in order to study also the temperature dependence of the critical velocity. Though 
we investigate the Green's functions it is important that they are coupled to the 
propagators describing collective motion and as a result their excitations hybridize 



3 Bogoliubov theory of a spin-1 Bose gas with a moving condensate 

The simplest possible approximation is the Bogoliubov approximation, which is valid 
at zero temperature or very close to it. In the Bogoliubov approximation the non- 
condensed atom fraction is assumed to be negligible, therefore the condensate atom 
number is approximately equal to the total number of atoms: Nc « N. 

There are two distinct phases distinguished by the number of condensate compo- 
nents, which we call as PI and P2 phases in analogy to the case of superfluid '^He. 
When the gas sample is maximally polarized all of the condensate atoms are in a 
single spin component and the condensate spinor is (f-^ — (1, 0, 0) |10lll| . In spite of 
the antiferromagnetic coupling the condensate is maximally magnetized due to the 
initial preparation and the conservation of the magnetization. When the polariza- 
tion of the gas is decreased the condensate wavefunction becomes multicomponent 
C^^ = Coj C-)- Note that the component Co = always in the Bogoliubov approx- 
imation, which is the consequence that the magnetization of the condensate points 
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Fig. 1. Dispersion relation of spin waves in the Bogoliubov approximation. The dashed (red) 
line is for vq = 0, the solid (green) line is for vo — «o,c ~ C2, the dotted (blue) line is for 
«o ~ 1-2 X vo.c- The parameters are T = 0, huJL = OAcsU, and Cs = 3 x 10~^c„. 



in the direction of the magnetic field (although our case it is just a Lagrange multi- 
plier). Consequently the x and y components of the condensate magnetization have 
to vanish and therefore Co = [TT. In this paper we restrict the discussion to the P2 
phase. 

The equation of state follows from the generalized Hugenholtz-Pines theorem [TD] . 
When the condensate is not fully polarized, that is when C_ ^ 0, we have two separate 
equations, one for the chemical potential and one for the magnetic field. 

M = eico + ncn, (8a) 
Hlul = mcs, (8b) 

where n = Uc — n.c,+ + n^^- is the total denstiy of atoms, and m = 71^^+ — n^^- is the 
total magnetization density. The condensate density in spin component r is given by 
nc,r = ncCr- Equations ([8| relate the chemical potential (fi), the magnetic field {tii^L) 
to the density n, the magnetization m and the velocity of the fluid (fiko/M). 
The self-energy in the Bogoliubov approximation is given by 

S^p(k,iuj„) = h^^ [{ncn + r m Cs)6ap&rs + ^M^7«Z7(c„ +r scs)] . (9) 

Combining Eqs. ([5]), ([6]), ([t]) and ^ and performing some straightforward algebra 
one obtains the following equations for the dispersion relation of the quasiparticles: 

= - p{cn + Cs)n + AcnCsUc^+nc- , (10a) 
p(^Ae + Ae^ = (uj - Ae^ (^uj + Ae^ . (10b) 

Where we have introduced the quantity p, which is provided by the solution of Eq. 



( 10a I. We have also introduced Ae = Ck — and Ae = e2ko-k — ekp- By measuring 



the wave number relative to the condensate k' = k — ko, the solution of Eq. ( 10b ) is 



cjk' = k' • Vo ± ^e^, +2pek', (Ha) 

from which we see that the speed of the quasiparticle for the condensate in rest, 
i.e. the proportionality constant in the starting linear part of the dispersion relation, 
is simply ci.2 — \J px^ijM . The quantity p is obtained from the quadratic equation 



(10a I 



Pi,2 ^ ± \^/{c^ - cs?n^ + 4e.nc,m2. (lib) 
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It is easy to see that both solutions (lib) are positive. The dispersion relation Eqs. 
([TT]) for Vq = are in complete agreement with the results of Ohmi and Machida [5]. 
The two solutions correspond to the density wave (with the + sign) and spin wave 
(with the - sign) excitations. When Cg ^ c„ the speed of the spin waves is proportional 
to (cs)^/^, while the speed of sound is proportional to {cnY^^- 



the first term in Eq. (11a 



The dispersion rel ation (11a) is very intuitive in accordance with Landau's idea: 
is just what we expect after a Galilean transformation to 
re condensate, while the second therm coincides with the 



the frame moving with t 
result in the frame where the condensate is at rest. 

When the condensate moves with velocity Vq, excitations with wavenumber par- 
allel to the condensate have a steeper dispersion while those which move against the 
condensate have a flatter one. The critical velocity vq.c is where the dispersion of the 
oppositely moving quasiparticles is horizontal. Since the spin wave velocity is smaller 
than the speed of sound, the instability due to the motion of the condensate appears 
in the spin wave dispersion first. Therefore the critical velocity is the speed of spin 
wave for a condensate at rest: vq^c = C2. In Fig. ^we show the spin wave dispersion 
for the three cases. 



4 Mean-field approximation at finite temperature 

In order to evaluate both the temperature and magnetization dependence of the crit- 
ical velocity we use the Landau criterion generalized to the spinor case according to 
the previous section. Here we use the finite temperature excitation spectrum of a 
spinor Bose gas in the P2 phase, worked out earlier in [TD], and evaluate the speed 
of spin wave of a condensate at rest, i.e. the critical velocity, at different temperature 
regions. At finite temperatures the condensate densities depend both on temperature 
and magnetic field. The partially polarized P2 phase occupies a finite region in the 
space of temperature and magnetic field (or rather magnetization) . When the mag- 
netization is zero the P2 phase goes directly to the noncondensed phase. However, 
when the gas sample is polarized the P2 phase first goes to a diefferent BEG phase, 
the PI phase, where the condensate spinor is single component. The phase diagram 
is depicted in Fig. [2] (a). This approximation takes into account the effect of thermal 
excitations in a dynamically consistent way: this approximation is conserving and 
gapless and also accounts for the hybridization of one-particle and collective excita- 
tions [9 . Note that to ensure the coincidence of one particle and collective motions 
already in the Hartree approximation one needs to sum up bubble diagrams 

For further discussions we consider the experimentally relevant situation, when 
Cs ^ c„. In fact, for ^"^Na the ratio of coupling constants is e = Cs/c„ « 3 x 10~^. Be- 
sides the thermal wavelength, we can introduce two other characteristic lengths: S.B = 
h/ {AMcsnc-Y^"^ is the mean-field coherence length, and ^' = MUbT / {ATrh^ric-) is 
the characteristic length governing the phase fluctuations of the order parameter. 
As discussed in [101, at sufficiently low temperature, where the longest characteristic 
length is the thermal wavelength A ^ C-BjC'> the excitation energies are linear in 
wavenumber: a;i^2 = \/ P\,ilMk' , where p\ = n{cn + Cs) — Acgn+n-ln is for con- 
ventional sound waves, and p2 = noCs -t- Acsfi +n^ / n describes spin waves. The zero 



temperature limit agrees with the the results (lib) and with [S], however here rij. is 
the total particle density in spin component r — {0, — } coming from the conden- 
sate contribution n^r and the thermal one nj,. Note that no = rip is entirely thermal. 
It is remarkable that at medium temperatures, characterized by ^j, <C vl, the square 
of the spin wave velocity becomes 

2 _ P2 _ 4nc, + ne,- Cs , , 
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Fig. 2. a) The phase diagram in the Bogoliubov-Hartree approxmiation. Bc{0) — Csn/{g^B) 
is the critical value of the magnetic field separating the PI and P2 phases at T = 0. b) the 
dependence of the square of the critical velocity (in units of CsTl/M) on the temperature and 
magnetic field in the P2 phase. 



but now an essential difference is that the densities appearing in Eq. ( 12 ) are those 
of the condensate, that is they exhibit significant temperature dependence. We do 
not discuss the critical region where the mean-field theory breaks down, anyhow, one 
expects that the critical velocity tends to zero there. Fig[2](b) shows the temperature 
and magnetic field dependence of the square of the critical velocity C2. The critical 
velocity vanishes at the P1-P2 phase boundary but even in the P2 phase its value can 
be very small when Uc,- is small. 

We have discussed the problem of critical flow in spin-1 Bose gases in the P2 
phase. We have shown that the critical velocity is the speed of spin waves and is 
much smaller than the critical velocity in scalar gases. We have not dicussed the PI 
phase, but one can expect, that there, close to the P1-P2 phase boundary the gap of 
the spin wave plays the role of critical velocity. 
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